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ABSTRACT 

An  experimental  investigation  was  performed  to  determine, 
in  situ,  the  change  in  mean  particle  size  across  the  exhaust 
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used  provided  excellent  results  within  the  exhaust.   However, 
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hampered  light  scattering  measurements  within  the.  motor. 
Particle  size  measurements  were  consistent  with  the  sizes 
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I.   INTRODUCTION 

The  importance  of  aluminum  oxide  (A£,?0-.)  particle  size 
produced  in  solid-propellant  rocket  motors  has  been  recog- 
nized since  aluminum  was  first  introduced  as  a  fuel  additive. 
Aluminum  added  to  a  propellant  will  increase  the  thrust 
produced  and  aid  in  suppressing  high-frequency  combustion 
instability.   Although  the  performance  improves  with  aluminum 
addition,  two  phase  flow  losses  adversely  affect  specific 
impulse  efficiency.   Specific  impulse  efficiency  is  the  product 
of  several  component  efficiencies,  i.e.,  chemical-kinetic, 
combustion,  two-phase  flow,  two-dimensional  flow,  submergence 
and  throat  erosion.   In  aluminized  propellants,  the  two-phase 
flow  loss  is  one  of  the  more  significant  losses.   Two-phase 
flow  losses  result  from  particle-gas  thermal  and  velocity  lags. 
The  latter  strongly  depends  on  particle  size.   Parametric 
studies  involving  small  motors,  with  throat  diameters  less 
than  one  inch,  at  AFRPL  [Ref.  1]  and  CSD  [Ref.  2]  have 
revealed  that  a  1  micron  change  in  particle  diameter  results 
in  a  change  of  1%  in  Isp.   The  effective  use  of  particle 
damping  for  suppression  of  combustion  pressure  oscillations 
is  also  strongly  dependent  upon  the  particle  size  distribution. 

One  of  several  analytical  computer  programs  used  to  pre- 
dict performance  of  a  solid-propellant  rocket  motor  is  the 
SPP  (Solid  Performance  Program)  [Refs.  3,4].   As  the  above 
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discussion  indicates,  it  is  of  critical  importance  to  accurately 
predict  particle  size.   This  is  especially  true  for  particle 
damping  with  minimum  smoke  propellants  which  typically  con- 
tain less  than  2%  metal  by  weight.   The  1975  version  of  the 
SPP  [Ref.  3]  contains  a  particle  sizing  empirical  model 
developed  by  Cohen.   This  model  was  based  on  correlating 
collected  exhaust  samples.   After  a  critical  review,  this  model 
was  considered  inadequate.   A  new  model  was  therefore  required 
for  a  new  SPP  version,  titled  Improved  SPP  [Ref.  4].   Inorder 
to  meet  SPP  accuracy  goals,  this  model  had  to  predict  particle 
size  within  ±10%.   Unfortunately,  there  are  no  adequate 
theoretical  models  relating  particle  size  to  propellant  and 
motor  parameters.   This  necessitates  the  use  of  empirical 
methods  for  correlating  particle  size. 

There  are  two  methods  that  can  be  employed  to  develop  an 
empirical  relation.   One  is  to  correlate  particle  size  data  to 
propellant  and  motor  variables,  as  was  done  by  Cohen.   The 
other  method  is  to  calculate  particle  size  based  on  a  critical. 
Weber  number  for  a  maximum  stable  drop  size  [Ref.  5].   After 
considering  the  relative  merits  of  each  method,  Hermsen  [Ref. 
5]  devised  an  empirical  model  based  on  collected  particle  size 
data  for  the  improved  SPP.   This  model  was  correlated  to  a  much 
broader  data  base  than  was  previously  available  to  Cohen. 

The  particle  diameter  correlated  was  D4 .. ,  the  mass-weighted 
diameter,  given  by: 
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d,  ,     =       y N.D4/  y N.D3 

4  3       L     l  1/    L     l  l 


Smith  [Ref.  6]  has  shown  that  D.^  was  the  most  appropriate 
average  to  use  in  predicting  two-phase  flow  losses. 

Hermsen's  empirical  relation  [Ref.  5]  is  presented 
below: 


D 


=   3.6404  D;2932(1_e-.00081635PT 


43      ------  -t 


where : 

D43  =  mass-weighted  average  diameter,  \i 

D  =  nozzle  throat  diameter,  in. 

^  =  ALO,  concentration  in  chamber,  g-mol/100  g 

P  =  chamber  pressure,  psia 

t  =  average  residence  time,  msec 

The  predicted  D.~>    is  now  applied  to  the  expression  below  in 
order  to  calculate  the  two  phase  flow  loss  [Ref.  7]. 


C  _  C . 

ETATP       =       C3UC4D    5/P*1:V08Dt6 


where 


ETATP   =   impulse  loss  effect,  % 

i  =      mole  fraction  of  condensed  species,  moles 

of  condensed  species  per  100  grams  of 
mixture 

D       =   D.-,  in  microns 
p  43 
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P     =   average  motor  chamber  pressure,  psia 

e      =   nozzle  expansion  ratio  at  ignition 
conditions 

D.     =   nozzle  throat  diameter  at  ignitiion, 
inches 

e     =   Naperian  base  2.71828 

C(  )   =   coefficients  that  are  a  function  of  D 

and  D 
P 

The  two-phase  flow  loss  is  applied  with  the  other  aforementioned 
efficiencies  to  a  theoretical  Isp.   The  Isp  loss  from  boundary 
layer  effects  is  then  subtracted  from  the  corrected  theoreti- 
cal Isp  to  obtain  delivered  Isp. 

Collecting  exhaust  products  is  feasible  only  for  small 
motors  and  the  techniques  used  often  result  in  variation  in 
the  measured  size.   Validation  of  the  SPP  model  must  be  based 
on  data  collected  within  the  motor  environment.   Light 
scattering  measurements  are  well-suited  for  this  application. 
These  measurements  are  more  easily  made  if  the  particles  are 
large  compared  with  the  wave  length  of  the  transmitted  light 
and  the  transmittance  is  greater  than  90%.   The  former  allows 
the  accurate  application  of  Fraunhofer  diffraction,  whereas  the 
latter  specification  is  necessary  in  order  to  satisfy  single 
scattering  requirements. 

The  Naval  Postgraduate  School  has  conducted  a  series  of 
investigations  to  determine  the  applicability  of  light  scatter- 
ing measurements  to  the  solid  rocket  motor  exhaust  products 
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[Refs.  8,9,10,11].  This  investigation  was  directed  at  modi- 
fying the  apparatus  and  techniques  developed  in  the  previous 
efforts  in  order  to  more  accurately  determine  particle  size. 
Validation  of  the  apparatus  was  also  a  goal.  The  SPP  particle 
prediction  model  can  then  be  verified  by  a  comparison  of  the 
predicted  size  to  the  experimental  data. 
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II.   THEORETICAL  BACKGROUND 

The  most  general  theory  of  light  scattering  was  developed 
by  Mie  and  is  presented  by  Van  de  Hulst  [Ref.  12].   Mie's 
theory,  although  not  restrictive  in  its  assumptions  in  regard 
to  particle  size  and  refractive  index,  is  however  very  com- 
plex to  utilize  for  data  reduction.   Dobbins,  et  al.  [Ref.  13] 
and  Powell,  et  al.  [Ref.  14]  derived  similar  expressions  for 
intensity  ratios  at  two  forward  scattering  angles  based  upon 
an  Upper  Limit  Distribution  Function  (ULDF) .  and  Fraunhofer 
diffraction.   The  ULDF  was  initially  proposed  by  Mugele  and 
Evans  in  Reference  15.   Although  the  expressions  developed  by 
Dobbins  and  Powell  are  extremely  useful  in  relating  intensity 
ratios  to  D-,2'  the  volume-surface  mean  diameter,  these  models 
are  still  somewhat  cumbersome  to  use  for  data  reduction 
purposes . 

Particle  sizes  found  in  solid-propellant  rocket  motors 
are  often  bi-modal .   Particles  in  the  larger  mode  are  the  most 
important  for  two-phase  flow  loss  calculations  and  are  gener- 
ally much  greater  in  diameter  than  the  wavelength  of  inci- 
dent light.   Scattering  by  these  large  particles  is  described 
by  Fraunhofer  diffraction.   Hodkinson  [Ref.  16]  and  others 
[Refs.  13,14]  have  demonstrated  that  diffraction  theory 
(within  the  forward  lobe)  can  be  used  to  approximate  the  Mie 
theory.   Buchele  in  Reference  17  summarized  the  experimental 
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techniques  for  determining  the  particle  size  by  measurements 
of  dif fractively  scattered  light.   Buchele  presents  a  function 
which  closely  approximates  the  curves  given  in  References  13 
and  14. 


1(6)   =   EXP [-(.57  a  8)2] 


where : 


1(0)   =   ratio  of  intensity  of  scattered  light  at 
some  angle  (0)  to  the  intensity  of 
scattered  light  at  theta  equal  to  zero 
degrees 

a     =   ttD/A  is  the  particle  size  parameter  for 
diameter  D-.2  an<3  wavelength  of  light 
lambda  ( A)  . 


The  above  expression  was  used  in  the  current  study  for  data 
reduction  to  obtain  D^2*   Unfortunately,  D.~  is  the  mean 
diameter  used  in  the  SPP.   Therefore,  D.,-  must  be  related  to 
D»~.    either  by  analyzing  exhuast  samples  or  by  approximate 
analytical  means. 

A  final  important  point  is  that  the  above  expression  is 
valid  for  a  <_  3  (the  center  lobe) .  For  the  expected  ranges 
of  theta  and  D->2'  this  proves  to  be  a  very  mild  restriction. 
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III.   EXPERIMENTAL  APPARATUS 

A.  ROCKET  MOTOR 

The  rocket  motor  used  in  this  investigation  was  the  same 
as  the  long  motor  described  in  Reference  11.   The  motor  com- 
ponents are  presented  in  a  photograph  in  Figure  3.1  and  the 
installed  motor  is  shown  in  Figure  3.2.   A  schematic  is  also 
presented  in  Figure  3.3.   The  propellant  grain  dimensions  re- 
mained the  same  as  in  the  previous  effort  [Ref.  11].   The 
grain  was  two  inches  in  diameter,  one  inch  in  length  and  had 
a  web  thickness  of  .725  inches.   The  grain  was  cylindrically 
perforated  and  allowed  to  burn  also  on  the  aft  end.   The 
igniter  was  composed  of  BKNO->.   Ignitiion  of  the  igniter  was 
achieved  by  heating  a  resistance  wire  with  a  12  VDC  source. 

Some  of  the  exhaust  products  were  collected  in  an  eight 
inch  diameter  stainless  steel  exhaust  tube.   The  products 
were  collected  from  both  the  side  walls  of  the  tube  and  the 
end  cap.   Another  particle  sample  was  collected  from  the  walls 
of  the  combustor,  near  the  nozzle  entrance.   Both  samples  were 
then  cleaned  and  mounted  for  observation  using  a  scanning 
electron  microscope. 

B.  LIGHT  SCATTERING  APPARATUS 

The  apparatus  shown  in  Figures  3.4  and  3.5  was  basically 
the  same  as  used  in  Reference  11  with  some  exceptions.   The 
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light  source  used  for  the  nozzle  exhaust  was  a  8-miliwatt, 
helium-neon  laser.   The  nozzle  entrance  also  used  a  8-milliwatt 
helium-neon  laser.   The  lasers  were  mounted  on  two  parallel 
optics  benches,  such  that  one  beam  could  pass  through  the  ex- 
haust and  the  other  beam  could  pass  through  the  motor  cavity 
adjacent  to  the  nozzle  entrance.   A  beam  expander/collimator 
was  used  with  the  helium-neon  laser  in  the  exhaust.   Each  beam 
was  intercepted  by  a  stop  placed  directly  in  front  of  the 
narrow  pass  filter.   The  narrow  pass  filter  was  used  to  eliminate 
extraneous  light.   The  stop  and  narrow  pass  filter  were  placed 
30.5  centimeters  from  the  motor  centerline  for  the  exhaust 
beam  and  11.5  centimeters  from  the  motor  centerline  for  the 
motor  cavity  beam.   The  exhaust  beam  narrow  pass  filter  had 
a  diameter  of  5.08  centimeters  and  the  motor  cavity  narrow 
pass  filter  had  a  diameter  of  2.3  centimeters.   Both  beams 
passed  through  the  centers  of  these  filters.   In  the  present 
arrangement  of  the  apparatus  this  resulted  in  a  maximum  angle 
for  collecting  scattered  light  of  4.8  degrees  in  the  exhaust 
and  6.0  degrees  in  the  motor  cavity.   The  placement  of  the 
optics  closer  to  the  motor  centerline  allowed  for  a  greater 
maximum  angle  to  be  recorded  than  in  the  original  design. 
Directly  following  the  narrow  pass  filters  were  condensing 
lenses.   Both  paths  used  a  lens  with  a  diameter  of  5.08 
centimeters.   These  lenses  had  a  focal  length  of  50  centimeters 
and  were  utilized  to  focus  scattered  light  onto  a  linear  array 
of  photodiodes.   The  array  consisted  of  1024  silicon  photodiodes 
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with  .025  millimeter  spacing.   The  apparatus  in  Reference  11 
positioned  the  lenses  such  that  only  angles  from  0  to  .04 
radians  could  be  observed.   In  actuality,  however  (as  a  conse- 
quence of  diffraction  around  the  beam  stop  and  beam  spreading) , 
the  range  of  useful  angles  where  only  scattered  light  occurred 
was  from  .02  radians  to  .04  radians.   In  the  current  arrange- 
ment the  diode  array  was  placed  1.0  centimeters  below  the  main 
beam  and  the  collecting  tube  diameter  was  enlarged.   This 
allowed  a  range  of  angles  from  .02  to  .07  radians  to  be 
observed. 

C.   DATA  ACQUISITION  AND  DATA  REDUCTION 

The  data  acquisition  system  is  presented  in  detail  in 
Reference  10.   The  system  controller  was  an  HP  9836S  computer 
and  the  HP  6942A  multiprogrammer  provided  A/D  conversion  and 
storage.   The  data  acquisition  system  was  programmed  to  take 
data  at  a  specified  point  in  the  motor  firing  by  inputing  a 
trigger  pressure  and  a  time  delay  after  the  trigger  pressure 
was  reached.   The  number  of  diode  scans  for  the  exhaust  and 
motor  cavity  were  eight  and  four  respectively  (limited  by 
the  current  data  acquisition  system  memory  size) . 

Two  methods  were  used  in  References  10  and  11  for  data 
reduction,  one  of  which  was  altered  for  this  study.   The  modi- 
fied method  pertained  to  the  iterative  graphics  technique. 
In  References  10  and  11  the  experimental  data  were  normalized 
and  compared  to  a  theoretical  profile  that  was  also  normalized 
to  some  assumed  scattered  light  intensity  at  zero  degrees 
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for  a  specified  D^2-   A  "best-fit"  of  the  experimental  profile 
to  the  theoretical  one  yielded  D-.~.   The  current  investigation 
determined'  that  this  method  resulted  in  ambiguity  in  deter- 
mining D^2  as  a  consequence  of  assuming  a  centerline  intensity 
for  the  scattering  profile.   The  method  used  in  the  present 
study  was  one  in  which  the  experimental  data  were  not  altered. 
A  theoretical  curve  was  generated  for  any  specified  value  of 
D^2  based  upon  a  selected  initial  angle  (9, )  and  the  correspond- 
ing intensity  at  that  angle  (I, )  from  the  unaltered  experimental 
data.   The  approximate  intensity  profile  for  a  polydispersion 
given  by  Reference  17  is: 


I2      =      I1EXP-{  (02  -  0^)  (.57    Dtt/A)2} 


Specification  of  a  series  of  values  for  0~  results  in 
values  for  I~  and,  therefore,  an  I  vs  0  theoretical  profile 
that  passes  through  the  experimental  data  at  I,  #  9]_.   If 
this  curve  does  not  fit  the  experimental  data,  then  other  values 
of  D-30  can  ke  chosen  and  a  new  theoretical  curve  generated. 
A  "best-fit"  method  is  subsequently  used  to  determine  D.,-. 
In  summary,  the  method  used  for  obtaining  the  mean  particle 
size  is  based  upon  fitting  an  approximate  (although  quite 
accurate  within  specific  angle  restrictions)  theoretical 
profile  for  a  polydispersion  to  the  unaltered  profile  obtained 
experimentally  over  a  small  range  of  scattering  angles.   The 
modified  data  reduction  program  is  presented  in  the  Appendix. 
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The  second  method  was  similar  to  the  first  and  also  used 
the  above  expression  solved  for  D^2 : 


D*2      =    -Ln(I2/I1)  (A/.57^)2/(92  -  6^) 


The  procedure  used  was  to  first  select  a  starting  angle,  0, , 
along  with  a  6-/6,  ratio.   A  value  of  I~  at  62  was  obtained 
from  the  experimental  data.   This  was  repeated  with  increasing 
values  of  92  and  then  again  with  different  values  of  6-/9-,  . 
The  results  were  plotted  graphically  as  particle  size  vs 
scattering  angle  for  each  angle  ratio.   A  nearly  constant 
value  of  D^2  vs  9  indicated  a  good  correlation. 
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TABLE  I 


T3 


K 


-n 


H 


Description 


Copper 
Nozzle 


Copper 
Nozzle 


Outside  Diameter 
(d,  inch) 


2.125 


2.125 


Length  (In,  inch) 


1.25 


1.25 


Throat  Diameter  (dt,  inch)      .28 


.23 


Slope  (converging) 
Angle  {y ,    degree) 


45 


45 
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Figure  3.1.   Motor  Components 
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TABLE  II 
LIGHT  SOURCE  SPECIFICATION 

A.  Helium-Neon  Laser  (Exhaust) 

1.  Manufacturer  Spectra-Physics 

2.  Model:  14  7 

3.  Type:  He-Ne  Class  IIIB 

4.  Output  Power:  8  mWatt 

5.  Beam  Diameter:  .92  mm 

6.  Beam  Divergence:  .87  nrad. 

B.  Helium-Neon  Laser  (Motor  Cavity) 

1.  Manufacturer:  Uni-Phase 

2.  Model:  1305P 

3.  Type:  He-Ne  Class  IIIB 

4.  Output  Power:  5  mWatt 

5.  Beam  Diameter:  .81  mm 

6.  Beam  Divergence:  1.00  mrad. 
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IV.   RESULT  AND  DISCUSSION 

A.  INTRODUCTION 

The  purpose  of  this  study  was  to  improve  and  validate 
the  accuracy  of  the  apparatus  used  for  particle  sizing  in 
the  rocket  motor  environment.   Also,  the  experimental  data 
will  provide  a  basis  to  verify  the  particle  size  prediction 
model  of  the  SPP. 

The  apparatus  had  been  calibrated  by  Harris  [Ref.  10] 
and  Kertadidjaja  [Ref.  11]  using  glass  and  polystyrene  spheres 
in  a  suspension  of  water.   Kertadidjaja  [Ref.  11]  also 
attempted  to  demonstrate  the  applicability  of  the  apparatus 
to  the  combustion  environment.   This  investigation  continued 
the  earlier  studies  with  several  refinements.   Three  propel- 
lant  types  were  used  in  this  investigation.   All  three  utilized 
a  GAP  binder  with  ammonium  perchorate  (AP)  oxidizer.   Each  had 
varied  aluminum  content,  with  the  percentages  present  being 
0%,  2%,  and  4.8%.   The  size  of  the  aluminum  powder  was  20 
microns.   The  major  ingredient  compositions  for  the  4.8% 
aluminized  propellant  were  14.7%  GAP,  45.7%  200  ym  AP,  24.6% 
25  urn  AP,  8.5%  TEGIN  and  4.8%  20  ym  A£  .   The  0%  and  2% 
aluminized  propellants  kept  the  same  solids  loading. 

B.  SYSTEM  CALIBRATION 

As  previously  mentioned,  the  apparatus  was  calibrated  by 
both  Harris  [Ref.  10]  and  Kertadidjaja  [Ref.  11].   The 
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apparatus  was  re-calibrated  in  this  study  as  a  result  of  the 
modifications  that  were  made.   The  standard  calibration  proce- 
dure described  in  detail  by  Harris  [Ref.  10]  remained  the 
same.   The  laser  beam  passed  through  a  glass  container  holding 
a  suspension  of  polystyrene  spheres  and  water  and  was  initially 
positioned  on  the  number  one  diode.   The  linear  diode  array 
was  then  moved  vertically  down,  with  the  aid  of  a  dial  indica- 
tor, such  that  the  number  one  diode  was  one  centimeter  below 
the  main  beam.   Since  the  diode  array  was  2.56  centimeters 
long  and  was  50  centimeters  from  the  focusing  lens,  this 
arrangement  allowed  for  a  range  of  measuring  angles  from  .02 
to  .07  radians.   Therefore,  the  scattered  light  striking  the 
diodes  would  not  be  affected  by  the  main  transmitted  beam  or 
the  light  diffracted  by  the  beam  stop  as  in  the  previous  two 
studies.   This  arrangement  was  used  both  for  the  exhaust  beam 
and  the  motor  cavity  beam. 

Two  polydispersion  samples  of  polystyrene  spheres,  with 
D,?'s  of  4.8  microns  and  10.2  microns,  were  used  for  the  cali- 
bration.  These  samples  were  chosen  because  the  D^.„  was  known 
and  the  expected  D-,2  °f  the  actual  aluminum  oxide  particles 
during  a  motor  firing  was  between  1  and  15  microns.   The 
results  are  presented  in  Figures  4.1,  4.2,  4.3  and  4.4.   These 
results  were  encouraging  in  regards  to  validating  the  modifi- 
cations and  the  method  for  data  reduction.   Calibrations  were 
also  conducted  within  the  motor  cavity  with  the  laser  beam 
and  scattered  light  passing  through  the  motor  windows.   Figures 
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4.1  and  4.2  are  representative  of  the  results  obtained.   This 
indicated  that  the  motor  windows  did  not  adversely  affect 
scattering  profile  and,  therefore,  bias  the  results. 

An  initial  point  of  concern  was  the  narrow  range  of  measur- 
ing angles  which  were  employed  in  the  present  apparatus. 
Since  large  particles  diffract  most  of  the  incident  light  at 
small  angles  and  very  small  particles  diffract  most  of  the 
incident  light  at  large  angles,  the  range  and  minimum  angle 
employed  could  not  have  resulted  in  sufficient  sensitivity  to 
particle  size  for  the  expected  range  of  particle  sizes.   Also 
shown  in  Figure  4.1  are  two  additional  theoretical  profiles 
for  D-,?  of  10  and  8.3  microns.   It  was  apparent  from  the  sensi- 
tivity of  the  slope  change  with  D^?  that  the  mean  particle 
size  could  be  determined  with  an  uncertainty  less  than  ±.5 
microns . 

The  "upsweeps"  in  the  curves  at  larger  angles  using  the 
two  angle  method  were  the  result  of  the  increase  in  slope  of 
the  scattering  profile  due  to  the  effect  of  truncation  by 
the  apparatus.   These  parts  of  the  curves  remain  only  to 
provide  a  reference  to  the  angle  ratios. 

Another  informative  point  in  apparatus  design  is  given 
by  Hirleman  [Ref.  18]  in  regards  to  placement  of  the  focusing 
lens.   Hirleman  notes  that  truncation  of  the  larger  scatter- 
ing anbles,  which  is  a  function  of  the  placement  of  focusing 
lens,  will  cause  a  biasing  against  the  small  particles.   The 
expression  given  in  this  discussion  is: 
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Z    =   (7Td/7.66A)D 
c 


where: 


Z    =   critical  distance  in  front  of  the  focusing 
lens  where  significant  energy  diffracted 
by  particles  less  than  diameter  D  is 
truncated  by  the  focusing  lens 

A    =   wavelength  of  the  incident  light 


For  a  given  Z  ,  all  particles  that  are . greater  than  D  in  size 
will  have  84  percent  of  their  diffracted  energy  pass  through 
the  focusing  lens.   In  the  exhaust  configuration  of  the  present 
apparatus  Z   was  equal  to  30.5  centimeters,  d  was  5.08  centi- 
meters and  A  was  .6328  microns.   This  yields  a  D  of  9  microns. 
Thus,  according  to  this  criterion,  some  of  the  scattered  light 
from  particles  below  9  microns  is  truncated  and,  therefore, 
results  in  a  bias  towards  larger  particles.   This  type  of 
biasing  was  small  in  the  present  apparatus  down  to  particle 
sizes  of  approximately  3  microns  (the  smallest  particles  in 
the  4.8  micron  D->~  polydispersion)  as  evidenced  by  the  good 
calibration  results  is  shown  in  Figure  4.3.   However,  the 
apparatus  could  be  improved  by  either  moving  the  lens 
closer  to  the  exhaust  (at  the  risk  of  apparatus  damage  from 
the  exhaust  jet)  or  by  using  a  larger  diameter  lens. 

Transmittance  effects  were  also  investigated  in  conjunc- 
tion with  the  calibration.   As  previously  mentioned,  trans- 
mittances  above  90%  are  required  to  meet  single  scattering 
requirements.   Calibrations  were  conducted  using  polystyrene 
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spheres  with  a  D^„  of  10.2  microns.   Scattering  data  were 
recorded  for  transmittances  of  85%,  70%,  60%,  50%,  and  30%. 
The  results  are  presented  in  Figures  4.5  through  4.14,  and 
summarized  in  Table  III.   Both  the  curve  fit  and  the  two- 
angle  methods  resulted  in  very  good  agreement  with  the  actual 
mean  particle  size  with  transmittances  as  low  as  50%.   Even  a 
transmittance  of  30%  did  not  severely  affect  the  measurement. 
These  results  indicate  that  particle  sizing  methods  based  upon 
single  scattering  theory  can  effectively  be  used  in  a  multiple 
scattering  environment  with  transmittances  significantly  less 
than  90%. 

A  test  was  then  conducted  to  determine  the  transmittance 
in  the  actual  motor  exhaust  using  the  2%  aluminized  propellant. 
The  procedure  was  the  same  as  the  one  used  in  Reference  11. 
The  transmitted  beam  was  directed  onto  the  diode  array  (with 
the  appropriate  filters  to  reduce  the  intensity)  and  the  beam 
profile  was  recorded.   The  motor  was  then  fired  and  the  profile 
was  again  recorded.   The  results  are  presented  in  Figure  4.15. 
The  results  indicated  a  transmittance  greater  than  90%,  with 
some  beam  spreading.   There  was  also  no  apparent  beam  shift. 
Therefore,  transmittance  loss  was  determined  not  to  be  a 
debilitating  factor  in  the  motor  exhaust  jet. 

The  calibration  curves  have  an  apparent  limit  of  approxi- 
mately .04  radians  for  which  scattered  light  can  be  recorded 
(for  example,  see  Figure  4.3).   However,  it  was  stated  that 
the  apparatus  was  designed  to  record  angles  up  to  .07  radians. 

33 


This  is  due  to  the  water-glass  interface  having  an  effective 
index  of  refraction  of  1.39.   Thus,  the  scattering  angles  are 
corrected  for  this  high  index  of  refraction,  resulting  in  a 
corrected  maximum  angle  of  approximately  .04  radians. 

C.   NOZZLE  EXHAUST  MEASUREMENTS 

One  test  was  conducted  with  a  nonmetallized  propellant. 
The  results,  presented  in  Figure  4.16,  show  no  scattering. 
Therefore,  the  scattering  observed  by  Kertadidjaja  [Ref .  11] 
at  smaller  angles  for  a  similar  test  was  probably  the  result 
of  beam  broadening  and  diffraction  around  the  beam  stop. 

Three  tests  were  conducted  with  both  the  4.8%  and  2% 
aluminized  propellants.   The  motor  firing  data  and  a  summary 
of  the  particle  size  measurements  are  presented,  respectively, 
in  Tables  IV  and  V.   From  light  scattering  measurements  each 
propellant  showed  a  similar  result,  one  test  with  a  D_~  of 
approximately  5  microns  and  one  with  a  larger  value  of  D^„ 
(7-12  microns).   This  test-to-test  variation  probably  resulted 
from  the  very  short  burn  times  and  the  corresponding  lack  of 
a  well-defined  steady-state  condition.   These  results  are 
presented  in  Figures  4.17  through  4.28.   Typical  S.E.M.  photo- 
graphs are  shown  in  Figure  4.29.   The  light  scattering  results 
were  quite  consistent  with  those  obtained  from  the  S.E.M. 
evaluation.   Although  the  sample  sizes  from  which  the  S.E.M. 
evaluations  were  based  were  relatively  small  (see  Table  VI) , 
it  still  provided  a  good  indication  as  to  the  validity  of 
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the  light  scattering  measurements.   It  should  be  noted  that 
D^2'  an<3  to  an  even  greater  extent  D.^,  is  strongly  dominated 
by  the  larger  particles.   One  or  two  large  particles  can 
significantly  change  D^2.   Future  S.E.M.  evaluations  should 
utilize  much  larger  sample  sizes. 

The  SPP  equation  for  D.-,  (which  is  primarily  a  function  of 
throat  diameter)  predicted  much  smaller  values  of  D.->  than 
observed  in  this  limited  set  of  data.   The  SPP  model,  however, 
is  based  on  particulate  data  from  motors  with  nozzle  diameters 
greater  than  1  inch. 

D.   NOZZLE  ENTRANCE  MEASUREMENTS 

Light  scattering  measurements  through  the  motor  at  the 
nozzle  entrance  were  unsuccessful  as  a  result  of  combustion 
light  existing  at  the  same  frequency  as  the  laser  light.   This 
was  the  case  for  both  the  helium-neon  laser  (A  =  .6328  ym) 
and  the  argon  laser  (A  =  .488  vim). 

The  analyses  of  the  S.E.M.  data  for  the  nozzle  entrance 
indicated  a  much  larger  D~2  and  D.^  as  compared  to  the  nozzle 
exit  measurements.   Therefore,  a  significant  reduction  in  size 
appears  to  have  occurred  to  the  larger  particles  across  the 
exhaust  nozzle. 
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TABLE  III 

EFFECT  OF  TRANSMITTANCE 

ON  MEASURED  PARTICLE  SIZE 

(Do2  from  Manufacturers  Data  Equals  10.2) 


Transmittance  Measured  D.,~ 

(%)  (microns) 


85  9.8 


70  9.6 


60  9.5 


50  9.4 


30  9.3 
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TABLE  IV 
SUMMARY  OF  MOTOR  FIRINGS 


Date  of   Wt  of    Press     Press     Burn    Residence    D+. 
Test    Alum      Pc     Max  P     Time      Time       (in.) 
(%)      (psig)    (psigX    (sec)      (msec) 


Oct  31,  85  0  350  366  3.8        -  .28 

Oct  22,  85  4.8  312  328  2.9  10.9  .28 

Oct  23,  85  4.8  328  359  2.3  11.3  .28 

Oct  18,  85  4.8  333  340  2.8  11.3  .28 

Nov  19,  85  4.8  450  466  .95  14.0  .25 

Oct  30,  85  2.0  296  296  2.8  10.4  .28 

Nov  26,  85  2.0  408  467  2.0  14.4  .28 

Nov  25,  85  2.0  566  575  1.8  16.7  .25 


Note:   Pressure  colume  labeled  P   indicates  pressure 
at  which  data  was  taken. 
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TABLE  V 
SUMMARY  OF  EXPERIMENTAL  RESULTS 


Wt  of    Press       D32         D32  D43         D43 

Alum      Pc     Light  Set.     S.E.M.       S.E.M.      SPP  Model 
(%)     (psig)     (psig)      Ext./Mtr.     Ext./Mtr. 


4.8       328       5.6       5.6/18.2      7.3/24.7      2.3 


4.8       333       5.5        -  /  -        -  /  -        2.3 


4.8      450      12.0      13.6/18.9     19.7/23.9      2.5 


2.0      296       7.0      12.2/14.4     18.6/19.6      2.3 


2.0      408       9.3        -  /  -        -  /  -        2.3 


2.0      566       4.6       4.9/10.8      6.0/13.2       2.5 


Note:   (-)  indicates  sample  was  not  taken  for  this  run. 
Diameters  are  in  microns. 
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TABLE   VI 
S.E.M.    PARTICLE    DISTRIBUTION 

Size  in  microns 
1       2       3       4       5     6     7       8     9     10     12     15     19     20     25     30 


Oct  23  Ext.  53     22     16     11       7     4      2        1 


Oct.    23  Mtr.  50     21     22     13       9     6     8       4     1       2       1       1  12 


Oct.   30  Ext.  45     42     11     13     21  2       2  1       1 


Oct.    30  Mtr.  24     15       3       6     28     3     7     13  6  3         12 


Nov.   19  Ext.  64     35     31       7       6     1     1       1  2       1  11 


Nov.   19  Mtr.  73     29       95     12     4111211         1  21 


Nov.  25  Ext.     44  39  22   3  14  4  1   1 


Nov.  25  Mtr.     37  30  11   1  15  8  3   8      3   4   3 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  study  have  shown  that  the  measure- 
ment of  D^2  in  rocket  motor  exhausts  can  be  made  using  quite 
simple  light  scattering  methods.   Apparatus  calibration 
results  and  S.E.M.  evaluations  of  the  collected  motor  exhaust 
products  support  this  conclusion. 

The  SPP  equation  model  underestimated  D.-.  for  this  inves- 
tigation.  This  indicated  that  the  model  predictions  cannot 
be  extrapolated  to  small  motors  with  throat  diameters  less 
than  one  inch.   The  light  scattering  method,  however,  can 
be  used  to  validate  the  SPP  model  for  the  larger  motors  upon 
which  it  was  based. 

In  order  to  improve  the  signal  to  noise  ratio  in  regards 
to  the  scattered  light  intensity,  it  is  recommended  that  the 
laser  beam  be  moved  closer  to  the  exit  plane  of  the  nozzle, 
where  the  particle  density  is  greater.   The  laser  beam  for 
this  study  was  7.3  centimeters  aft  of  the  nozzle  exit.   To 
resolve  the  combustion  light  problem  it  may  be  necessary  to  use 
a  different  laser,  such  as  a  C0~  that  has  a  different  wavelength 
than  present  in  the  combustion  process.   Also,  other  binder 
systems  should  be  tested  to  determine  if  they  significantly 
alter  the  light  spectrum. 
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APPENDIX 
PROGRAM  LISTINGS 


10  •  ** *******  ***x****xx»»******************»****»*****x*xx**x*xxx*xx*x* 

20  ixx**x»mx*x****xx*         RDC  ***********xx***x*!(* 

30  i***x*****x**xxx**«*      PLOTS  RAW  DATA      ***xxx***xxxx******* 

40  i******x**xix*<xx***       AVERAGES        ix**»*xx***xx**«xxx» 

50  !******x*****x*x****       FILTERS        **xx*xx****x*x*x***» 

60  i*******************  DETERMINES  HFAN  DIAMETER  *xx*»*xx*x*x ******** 

70  !****x»*»*x***x*xx*»  BY  INTERACTIVE  GRAPHICS   ******************** 

80  !*<*****************    AND  THE  TWO-ANGLE  METHOD    *xx*x*xx***xxx**x*** 

90  ! »*x**#***»**xx»**»*   Robert  Kelly  Harris  84  ******************** 

100  i***»*****»**»xxx**»   Revised,  John  S  Rosa  fl5  *x*x*ixx****x***x»xx 

110  i ******** x***x»******************************»****************»****x 

130  OPTION  PA^E  1 

130  TOM  /Hrdqaus/  Av2(1024) 

131  COM  A2M024)  PUFFER 
140  COM  /Gauss/  TH1024),L 

150  CC^  /Max/  M?,M5,Xt,Yt;Xw,Y«,Xx 

160  COM  /Readata/  B,P,H,Q3$l201,B4$[20],Zz$I2B],Yl(B192)  BUFFER  ,Y2(8192)  BI'.FFE 

R 

170  COM  /Two/  AvlU024),M,Ml,F 

18?  DIH  Scans(R) ,XU0?4) 

190  COM  /Plots/  Pl<(201,P2$[2(ll,P3$[20],P4i[2fl] 

200  INTFGER  Qraf (1:12480)  BUFFER 

?10  Choose:  PRINT  CHRSU2) 

220  Old=0 

230  PRINT  "  TO  LOOK  AT  ANY  DATA  FILE  THE  PROGRAM  NEEDS  SOME  STARTING  INFORMATI 

ON" 

240  PRINT  " 

250  P'WT  "       TO  ACCOUNT  FOR  THE  CHANGE  IN  WAVELENGTH  IN  THE  MEDIUM" 

260  PRINT  "" 

270  PRINT  ■  ENTER  THE  INDEX  OF  REFRACTION  OF  THE  MEDIUM" 

280  PPINT  "" 

2?8-  PRINT  "  WATER  =  1.33" 

300  PPINT  "  AIR  =  1.0" 

310  PRINT  "  ESTIMATE  OF  EXHAUST  =1.1" 

320  PEEP 

330  INPUT  "       THIS  VALUE  ADJUSTS  THE  COMPARISON  CURVE  TO  THE  MEDIUM" ,M1 

340  PRINT  CHR$(12) 

350  PRINT  ■   TO  ACCOUNT  FOR  REFRACTION  OF  LIGHT  AT  BOUNDARIES  BETWEEN" 

360  PRINT  »■ 

370  PPINT  '  THE  MEDIUM  AND  AIR  YOU  MUSI  ENTER  THE  INDEX  OF  REFRACTION   " 

380  PRINT  " 
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390  PRINT  "     OF  THE  COMBINATION  OF  THE  MEDIUM  AND  ITS  BOUNDARY" 

400  PRTNT  ■■ 

410  PRINT  "       THIS  VALUE  IS  APPLIED  DIRECTLY  TO  THE  DATA1 

430  PRINT  " 

430  PRINT  "         AIR=1.0  PLEXIGLASS  *  MATER  =  1.39" 

440  PRTNT  " 

450  PRINT  'ESTIMATE  OF  MOTOR  CAVITY  COMBUSTION  PRODUCTS  i  UINDOU  =  1.22' 

460  PRINT  " 

470  PRINT  '  ESTIMATE  OF  ROCKET  EXHAUST  =  1,1  OR  1.0" 

480  INPUT  H 

490  PRINT  CHR$(12) 

500  INPUT  ■  ENTER  LASER  WAVELENGTH  (HeNp=.6328,Ar=.488)*,L 

510  INPUT  'ENTER  THE  DESIRED  PLOTTING  INTERVAL  OF  DIO^E  ARRAY  (2,4.6)*,H 

520  F=?00  !nn  FOCUL  LENGTH  OF  OBJECTIVE  LFNS 

521  XU)=1 

522  Tl(1)=,02 

530  Diod=J25  Inn    DIODE  SPACING  MAY  BF  ,03 

540  FOR  1=1  TO  1023 

550      X(I+1)=I+1  ICREATF  A'J  ARRAY  OF  DTODF  NUMBERS 

560      THI+1)=ATN((I»Diod)/(F«M))4Tl(l)      "CREATE  AN  ARRAY  OF  THETA 

570      NEXT  I 

580  PRTNT  CHR$(12) 

590  CALL  DisplavKOld)  « JUST  PRINTS  INITIAL  REMARKS  ON  CRT 

600  IF  Old  THEN  CALL  Rpvipw(AvK*))      'THESE  TWO  LINES  ALLOW  FOR 

610  IF  Old  THEN  GOTO  Gauss  (REVIEWING  REDUCED  DATA  FILES 

620  PEEP 

630  INCUT  '         INPUT  TWO  FILENAMES  NOW  (NO-PART  ,PART)\Q3$,04* 

640  INPUT  "      ENTER  '1'  FOR  8  SCANS    '2'  FOR  4  SCANS\P 

650  7;i=MNlFRNAL,4,r  !  USES  LEFT  DISK  DRIVE  FOR  RAW  DATA  FILES 

660  CALL  Rpadata 

670  CALL  Shift 

680  BEEP 

670  Serpen:  ! 

"08  PPTNT  CHR$(12)  ICLEAR5  SCREEN 

710  Pi$="VOLrAGE" 

720  P2t="DinDF  NUMBER"  !  STRINGS  FOR  PLOTS 

730  P3$-"VnLTAGE  ws,  DIODE" 

740  P4$="RAW    DATA" 

750  M5=1024 

751  Mfl=0, 

760  Xt=64  !   SET  UP  VALUES  FOR  PLOTS 

770  Yt=,04 

780  Xn=4 

790  Y*=5 

791  Xx=0 

800  CALL  Plot(Y2(»))  !  DRAWS  AXES  ,FTC... 

810  GSTORE  Graft*)  !  STORES  GRAPHICS  DT SLAY  JUST  HADE 

820  r/LL  Dataplot(B,YK»),H)  IPLOT  NO-PARTICLES  DATA 

830  PRINT  "  NO-PARTICLES  PLOTTED" 

840  PRINT  " 

850  PRINT  '  HOW  MANY  SCANS  SEEM  TO  BE  GOOD  7" 

860  INPUT  J 
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0  PRINT  " 

1  PRINT  '        WHICH  SCANS  ARE  GOOD  ?...ie.  ..1,2, 4,5, 7, INCLUDE  LAST  CO 

HMA"  "LAST  COMMA  IS  REQUIRED  OR  YOU  HAVE  TO  HIT  CONT..  TWICE 

890  INPUT  ScansU) 

9J0  PRINT  CHRM12) 

919  PRINT  USING  ■/////////' 

920  PRINT  '  AVERAGING  THE  SELECTED  SCANS" 
938  GRAPHICS  OFF 

948  CALL  A»erage(J,Scans<*),YK*),Av1(i)) 

950  PRINT  USING  ,B" 

960  GLOAD  Graf(i)     ILOADS  GRAPHICS  ARRAY  RATHER  THAN  WASTE  TIME  RE-DRAUING 

978  GRAPHICS  ON 

988  CALL  Result(AvK*),X(*),H)  !  PLOTS  1024  ELEMENT  ARRAYS 

990  PRINT  USINC  ■//' 

1080  PRINT  ■  Average  Intensity  No-Particles" 

1010  ON  KEY  0  LAPEL  "AVERAGES  GOTO  Serpen 

1028  ON  KEY  1  LABEL  "PLOT-PARTICLES"  GOTO  1060 

1830  PRINT  "  PRESS  KEY  *  8  TO  RE-AVERAGE  OR  *  1  TO  CONTINUE" 

1040  Standby:  !  MANUAL  CALLS  THIS  INTERRUPT  DRIVEN  PROGRAMMING 

1050     GOTO  Standby  !  LOOPS,  WAITING  FOR  USER  TO  DECIDE 

1060  PRINT  USING  "?" 

1070  OFF  KEY  1  ! HELPS  AVOID  CONFUSION  BY  CLEARING  THAT  POX 

1080  CALL  Dataplot(B,Y2U),H)  ! PLOT  PARTICLES  DATA 

1098  PRINT  USING  "/" 

1100  PRINT"  PARTICLE  DATA  PLOTTED" 

1110  PRINT  "       FOR  A  HARD  COPY  OF  THIS  RAW  DATA  PRESS  KEY  ♦  6" 

1128  PRINT  "  " 

1130  PRINT  "  OR   TO  CONTINUE  PRESS  KEY  *  1" 

1148  ON  KEY  6  LABEL  'HARD  4  RAW"  GOTO  Raw 

1158  ON  KEY  1  LABEL  "CONTINUE  "  GOTO  Select 

1160  GOTO  Standby 

1178  Raw:  CALL  Plot(Y2(»),l)  !  THE  ONE  (1)  IS  AN  OPTIONAL  PARAMETER 

1180     CALL  Dataplot<B,Y2(*),H)     !  WHICH  IS  USED  TO  GET  HARD  COPIES 

1190  Select:! 
1200     OFF  KEY  6 
1210     OFF  KFY  1 

1220  PRINT  CHR$(12) 

1230  PRINT  USING  '////" 

1240  PRINT  "  HOU  MANY  SCANS  SFEM  TO  BE  GOOD  ?" 

1250  INPUT  J 

1260  PRINT  "  WHICH  SCANS  ARE  GOOD  ?.. ,ie... 1,2,4, 5, 7, INCLUDE 

LAST  COMMA" 

1270  INPUT  Srans(») 

1280  PRINT  CHR1U2) 

1290  PRINT  USING  "/////////■ 

1308  GRAPHICS  OFF 

1310  PRINT  "  AVERAGING  THE  SELECTED  SCANS' 

1320  CALL  AverageU, ScansU), Y2(»),Av2(«)) 

1330  PRINT  USING  "0" 

1340  P4$="AVERAGED  SCANS' 
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1350 
1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430  PRINT 

CONTINUF" 


CALL  Plot(Av2(»)) 

GRAPHICS  ON 

GSTORE  Graf(») 

CALL  Result(Av2(«),X(*),H) 

PRINT  USING  7/" 

PRINT  ' 

ON  KEY  2  LABFL  "A'JFRAGES"  GOTO  10B0 

ON  KEY  3  LABEL  'SUBTRACT'  GOTO  1450 


Average  Scattered  Intensity" 


HIT  KEY  *  2  TO  RE-AVERAGE  OR  KEY  ♦  3  TO 


1440 
1450 

1460 

1470 

1430 

1490 

1500 

1510 

o-PartulPSH 

15?0 

1530 

1540 

1550 


GOTO  Standby 

PRINT  USING  'P" 

OFF  KEY  3 

MAT  Avl=  Av2-Avl 

Gl.OAD  Craf<») 

CALL  Resolt(Avl(*),X<*),H) 

PRINT  USING  "///" 

PRINT  • 


"SUBTRACTS  NO-PARTICLES  FROM  PARTICLES 


Plot  of  the  Difference  Between  Particles  and  N 


PRINT  ■ 
PRINT  ' 
ON  KEY  6  LABEL 


KEY  *  6  FOR  HARD  COPY' 
KEY  *1  10  CONTINUE  " 


ON  KEY  1  LABEL 

1560  GOTO  Standbv 

1570  CALL  Pl«t(Avl(*),l) 

1580  CALL  Result(Avl(*),X(*),H) 

1590  PLOTTER  IS  3, 'INTERNAL' 

1600  PRINT  USING  '//' 

1610  PRINT  ' 

1621  PRINT  " 

1630  PRINT  ' 

1640  PRINT  ' 

1650  INPUT  ' 

1660  IF  Fil=0  THEN  Gau?? 

1670  P4t='FILTEREI>  DATA' 

1680  CALL  Filter(Avl(*),Fil) 

1690  PRINT  CURt(12) 

1700  CALL  Plot(Avim) 

1710  IF  Old  THEN  MAT  X=  Tt 

1720  CALL  Result(Ai>l(*),X(*),H) 

1730  PRINT  ' 
o-Partulp?" 

1740  PRINT  ■ 
TEP" 
1750 
1760 
1770 
1790 
1790 
1800 
1810 


HARD  AVERAGE"  GOTO  1570 
FILTER  ■  GOTO  1600 


ENTER  THE  NUMBFR  OF  TIMES  YOU  WISH  TO  APPLY' 
THE  DIGHAL  FILTER  FOR  SMOOTHING" 
EXAMPLE********  10" 
TAKES  ABOUT  1.5  MINUTES" 
ZERO  IS  ALSO  ACCEPTABLE", Fil 


Plot  of  the  Difference  Between  Particles  and  N 


PRINT  " 

PRINT  " 

PRINT  ■ 

PRINT  ' 

PRINT  ' 

ON  KEY  8  LAPEL  "HARD  FILTERED' 

ON  KEY  9  LABEL  '  NORMALIZE  ' 


AFTER  APPLICATION  OF  A  DIGITAL  FIL 

for  a  HARD  COPY  ■ 
PREPARE  the  PLOTTER  and* 
PRFSS  KFY  t  9" 
OR" 
PRESS  KEY  »  9  TO  NORMALIZE" 


GOTO  1838 
GOTO  Gauss 
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1820  GOTO  Standby 

1830  CALL  Plot(AvK»),l)                     !  The  1  is  for  hard  copy 

1840  CALL  Res»lt<Av1<*),X(*),H> 

1850  OFF  KEY  9 

1860  Gauss:  OFF  KEY  2                     !KEY  1  OPTION  IS  STILL  VALID 

1870  OFF  KEY  8 

1880  ON  KEY  1  LABEL  '   FILTER  ■  GOTO  1600  !  ALLOWS  ONE  TO  SMOOTH  OLD  DATA 

1870  !  WHEN  REVIEWING  HEWED  FILES 

1891  01d=l   !  THIS  ALLOWS  FILTERING  EVEN  AFTER  ENTERING  GRAPHICS  MODE.  THE 

1892  i  DIODE  COORDINATES  AREN'T  NEEDED  ANYMORE  *SFE  LINES  NEAR  CALL  FILTER 
1900  PRINT  CWW12) 

1910  'Av2  Array  is  Normalized  PUT    AvlU)  is  Still  Saved  For  Re-work  if  Needed 

1920  MAT  Av2=  AvI       •  TO  BEGIN  WITH  THE  ARRAY  IS  ASSUMED  TO  BE  NORMALIZED 

1930  i  AND  THE  USER  CHANGES  THIS  WITH  INTERACTIVE  GRAPHICS 

1940  i  IN  THE  SUBROUTINE  "Compare" 

1950  P1«="FILTERED  VOLTAGE 

1960  P2$='THETA  (rad)"                                                               (STRINGS  FOR  PLOTS 

1970  P?*="IN!ENSITY  vs.THETA" 

1980  P4t="CURVE  FIT  RESULTS' 

1990  M5=,070                                                                                 !SET  UP  VALUES 

2000  Xt=.005                                                                                "FOR  PLOTS 

2010  Yt=,l 

2920  X»=2 

2030  Y«=2 

'2031  Xx=.n2 

2040  CALL  Plot(AvH»5) 

2050  CALL  Res!ilt(Avl(«),TK»),H)               lpl-0T  0F  NORMALIZED  INTENSITY  vs.  THETA 

2060  GSTORE  Graff*)                                     !  SAVES  THE  SCREEN  IMAGE 

2070  CALL  CoHp.ve(A»l(*),M,Hl>Graf(*>> 

2080  GRAPHICS  OFF 

2098  ON  KEY  4  LAPEL  "OTHER    ARRAY'  GOTO  Choose 

2100  ON  KEY  5  LABEL  '  TWO-ANCLE"  GOTO  Bvchele 

2120  ON  KEY  6  LABEL  "STORE  DATA"  GOTO  Storedata 

2130  OH  KEY  7  LAPEL  '      QUIT      '  GOTO  Quit 

2140  CALL  Hencel                                                        !  PRINTS  OPTIONS  ON  THE  SCREEN 

2150  GOTO  Standby 

2160  Buchele:     " 

2170  CALL  Two^gle 

2180  CALL  Merrnel 

2190  GOTO  Standby 

2200  Storedata:' 

2210  CALL  Store 

2220  CALL  Mennel 

2230  GOTO  Standby 

2240  Quit:     ! 

2250  END 

2260  SUB  Averane(J,Scans(*),Y(*),Av(*))     ! AVERAGES  SELECTED  SCANS 

2270  MAT  ftv=  (0)                        'INITIALIZES  THE  ARRAY  LOCAL  TO  THIS  ROUTINE 

2280  FOR  1=1  TO  J  STEP  1 

2290  K=(Scans(I)-l 1*1024+1                   !  THIS  COUNTER  IS  THE  BEGINNING 

2300  L=Scans(I)»1024                            !  AND  THIS  ONE  THE  END  OF  BLOCKS 

2310  FOR  I1=K  TO  L                       'OF  1024  IN  THE  OVERALL  DATA 
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2320 
2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

24?0 

2480 

2490 

2500 


2530 
2540 
2511 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2631 
2640 
2650 
2660 
2670 
2680 
2700 
2710 
2720 
2730 
2740 
275B 
2760 
2770 
2790 
2790 
2800 
2810 


I2=U-K+1 

Av(I2)=Av<I2)+Y(I1) 

NEXT  11 
NEXT  I 

MAT  Av=  Au/(J) 

SUPFND 
SUP  Dataolot(B,Y(*),H) 
COM  /Max/  M7,M5,Xt,Yt,X«,Y«,Xx 
IDIR  0 
LORG  4 

FOR  J=l  TO  B 

MOVE  5e0+J*5O,M7-,05 


\»m 


NTFR  IS  ALWAYS  BETWEEN 


!B  IS  4  OR  8 (THE  NUMBER  OF  SCANS) 


LAPEL  J; 
LINE  TYPE  1 
K=(J-1)M02441 
L=J*lfl24 
MOVE  1,Y<K) 

FOR  T=K  TO  L  STEP  H 
M-K+1 
PLOT  X, Y( 1) 
NEXT  I 
NEXT  J 
PENUP 
M8=0 
SUBEND 

SLIP  Res«]t(Iaw(*),D(»),H) 

CLIP  ON 

MOVE  DU),Iav(1) 

FOR  1=1  TO  1024  STEP  \ 
DRAW  D(I),Iav(I) 
NEXT  I 
PENUP 
SUBEND 
Spit:     ! 

SUB  Plot (Y(*), OPTIONAL  Hard) 
COM  /Max/  M7,M5,Xt,Yt,X*,Y",Xx 


'EACH  SCAN 
!   SEE  NEXT  LINE 
!  HELPS  KEEP  TRACK  OF  FACH  SCAN  AS  IT  APPFARS 

! BEGINNING  OF  EACH  SCAN  AND 
"THE  END  WITHIN  THE  TOTAL  BLOCK 
!  MOVE  TO  THE  FIRST  POINT 

!  THIS  GIVES  1  TO  11124  FOR  ABSCISSA 


!THE  AVERAGE  INTENSITY  IS  PLOTTED 


i IF  THE  OPTIONAKHardHS  PFCE1VFD 
i THE  FIGURE  GOES  10  THE  PLOTTER 


COM  /Plots/  P1$I201,P2$[201,P3$120],P4M20] 


M7=MAX(Y(»))»1 

M6=.1*H7 

G1NIT 


SELECT  NPAR 

CASF  1 

PLOTTER  IS  3,'INTERNAL' 

CASE  2 

PLOTTER  IS  705/HPGL1 
END  SELFCT 
GRAPHICS  ON 
CSIZE  5, .6 
DEG 
LDIR  0 


!    A  SCAI  FNG  VARIABIE 
!  ANOTHER  SCALING  VARIABLE 
'G1NIT  IS  JUST  GOOD  PRACTICE  SO  YOU  KNHW 
'WHERE  YOU  ARE  BEGINNING 
'THIS  DFTFCTS  IF  IHE  HARD  TOPY  IS  DESIRED 
!IT  COULD  BE  DONF  WITH  IF  THEN  LOGIC 
!BUT  IS  PRESEN1ED  FOR  FAMILIARIZATION 
ISINCE  IT  IS  MORE  POWERFUL  IN  COMPLEX 
"SITUATIONS 

!  TO  BE  ABLE  TO  SEE  THE  PLOT 
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2820  LO»G  5 
2B30  MOVE  75,9*5 

2840  LABEL  USING  'K";P4« 

2850  MOVE  75,91 

2860  LAPEL  P3* 

2870  HOVE  75,20 

2880  LABEL  P2$ 

2890  LDIR  90 

2990  MOVE  18,60 

2918  LABEL  PI* 

2920  VIEWPORT  30,125,30,35  !  SCREEN  UNITS  FOR  MARGINS 

2930  FRAME  !  DRAMS  A  BOX 

2940  WINDOW  Xx,H5,-M6,M7 

2950  AXES  Xt,M6,Xx,0,Xii,l,5 

2960  LDIR  0 

2970  CSIZE  3,. 5 

2980  LORG  8 

2990  CLIP  OFF 

3000     FOR  I=-M6  TO  H7  STEP  M7/10  "NUMBER  THF  Y  AXIS 

3019  MPVE  Xx,I 

3020  LAPEL  USING  '# ,MD . DDDD" j I 

3039  NEXT  I 

3040  CSIZE  3M6 
3050  LDIR  90 
3060     LORG  8 

3070     FOP  I=Xx  TO  H5  STEP  Xt  iNUMBER  THE  X  AXIS 

3080     MOVE  I, -MA 

3090     LAPEL  USING  "t,K"}I 

3100     NEXT  I 

3110     PENUP 

3120  SUBEND 

3130    SUB  Co«p^rp(Avl(»),H,M, Ml, INTEGER  Graft*)) 

3140    COH  /Ga»ss/  T1(*),L 

3150    COM  /Hrdqaus/  Av2<») 

3160    COH  /Max/  M7,M5,Xt,Yt,X«,Yn,Xx 

3201  Poot:  ! 

3202  PRINT  ■  INPUT  REFERENCE  ANGLE  ■ 
3293    INPUT  Tt 

3295  FOR  N=l  TO  1024  STFP  1 

3296  IF  T1(N)>=Tt  THEN  3208 
3207    NEXT  N 

32P8  lt=0 

3209  FOR  I=N-5  TO  N+5 

3?10  It=It+AMl(I) 

3211  NEXT  I 

3212  Al=It/ll 

3220       OFF  KEY  "GETS  RID  OF  ALL  LABELS  ON  KEYS 

3239  ON  KEY  3  LABEL  "MENUE"  GOTO  4160 

3240  PLOTTER  IS  3, 'INTERNAL'      !  IN  CASE  A  HARD  COPY  WAS  JUST  MADE 
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3250 
3260 

Hard=0 
Cent<?rline=l 

!  (0)  SO  ONF  DOESN'T  EXIT  TOO  SPUN 
•  THE  INITIAL  NORMftllZING  VALUE 

3270 

D=10. 

!  INITIAL  PATICLE  MEAN  DIAMETER  IN  MICRONS 

3280  Change:    CSIZE  4,. 6 

3290 

PRINT  USING  ■//////' 

3300 

PRINT  ' 

STAND  BY  FOR  CURVE" 

3310 

HH'v'E  0,1 

3320 

CLIP  ON 

3331 

DIM  T(l:1024) 

3340 

FOP  1=1  TO  1024  STEP  H 

1 

3350 

T<I)=Al*EXP(-DA2*((TKI)A2-TtA2)»(((Ml».57*PI)/L)A2))) 

3370 

DRAW  T1(I),T(I) 

3375 

IF  <T1(I)*((PI*D*I11)/L))>3.8  THEN  3440 

3380 

NEXT  I 

3440 

PENDP 

3450 

IF  Hard  THEN  3550 

•EXITS  ROUTINE  IF  A  HARD  COPY  WAS  JUST  MADE 

3469 

PRINT  USING  "/////////" 

3470 

PRINT  ■ 

USE  THE  KNOB  TO  VARY  THE  PARTICLE 

SIZE' 

3480 

PRINT  " 

OR  HIT  KEY  *  6  FOR  HARD  COPY" 

3490 

PRINT  ■  " 

3500 

PRINT  " 

KEY  *  9  ALLOWS  YOU  TO" 

3510 

PRINT  " 

SELECT  NEW  REF  ANGLE" 

3520 

PRINT  ■  ■ 

3530 

PRINT  ■ 

HIT  KEY  *  3  TO  GET  OUT" 

3540 

PPIMT  USING  ■//////" 

3541 

OFF  KEY  0 

3550 

ON  KEY  6  I  ABEL  'HARD  COPY" 

GOTO  Hardgauss 

3560 

ON  KEY  9  LAP^L  "NORMALIZE" 

GOTO  Poot 

3570 

ON  KNOB  .05  GOTO  Pulse 

!(.B5)   IS  TIME  INTERVAL  IN  WHICH 

3580  Haiti    GOTO  Wait 

!  PUI  SES  FROM  THE  KNOB  APE 

3598  Pt 

ilse:  PRINT  USING  "g" 

!  COUNTFD  AND  THIS  NUMBER  IS 

3600 

Strnqll="Hean  Diameter  ' 

!  USPD  PY  THE  INTERRACTIVE 

3610 

Strng2*="  Microns" 

!   GRAPHICS  TO  VARY  THE  PARTICLE 

3620 

Count=KNPPX 

I  SIZE  AND  PLOT   THE  ASSOCIATED 

3630 

D=DR0UND(D+Count/15,2) 

!   GAUSSIAN  APPROXIMATIONS  OF 

3640 

GLOAD  Graf(») 

!  SCATTERING  PROFILES 

3650 

LDIR  0 

3660 

LORG  8 

3670 

MOVE  M5,.8*M7 

3680 

LAPEL  SfrngmVALKDHStrni 

|2* 

3690 

GOTO  Change 

3701  Hardqauss: ! 

3710 

PRINT  " 

PREPARE  The  PLOTTER- 

3720 

PRINT  ■ 

PRESS  CONTINUE  for" 

3730 

PRINT  ■ 

HARD  COPY" 

3740 

PAUSE 

3750 

Hard=l 

!  SO  THAT  SUBEXIT  OCCURS  AFTER  HARD  COPY 

3768 

CALL  Plot(Av2(*),l) 

!   liftHARD  COPY 
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377(1  CALL  Result(Av2<»),T1(*),H) 

3780  PRINT  USING  '?" 

3790  LORG  8 

330  0  MOVE  M5,.8»M7 

3810  LDIR  0 

3820  LABEL  Strpql«VAL$<D)4Strnq2i 

3840  GOTO  Change 

4160  Subexit :    ! 

4170  SUBEND 

4180  SUB  Menuel 

4190  PRINT  USING  '8' 

4200  PRIN1  "                YOU  CAN  RF-HQRMALIZE" 

4210  PRINT  "*#x»m#MK#XM»x»***»mm**m**mH*x*mxm^ym**»xx1' 

4220  PRINT  "                 YOU  CAN  RE-AVERAGE  (New  Data  (fob)' 

4230  PRINT  "»*»*MM*«MM»**mm*M#*x»*xxxiHfxmx#mxx**»x****j{»xxn 

4240  PRINT  '              PRESS  KEY  *  4  TO  LOOK  AT  OTHER  DATA' 

4250  PRINT  "                     OLD  OR  NHJ  ■ 

4260  PRINT  '                   MOTOR  /  EXHAUST* 

42"8  PRINT  "*#»»x»»x»«»»k«#xxx»«xx*xx»***»xx* *mx*MXXM*xit»»xx *m**x »" 

4280  PRINT  "                                               PRESS  KEY  *  5  ' 

4290  FRINT  ■                                     FOR  THE  '   TWO-ANGLE   '  METHOD" 

4300  PRINT  '                                          OF  PARTICLE  SIZING" 

4310  PRINT  "***»m»»»m*»mm¥<mmKxmmx***y»mm**xmx*ifxx-" 

4320  PRINT  "          TO  STORE  THE  REDUCED  DATA  PRESS  KEY  *  6 

4338  PRINT  "m»»#*«x*«»*»»»xmMmm»*x**xmxxmxxrxx«mx*mx*xx" 

4340  PRINT  "            TO  END  THIS  SESSION  PRESS  KEY  *  7" 

4350  SUPEND 

4360  SUB  DisplavKOld) 

4370  PRINT  "                 TO  REDUCE  NEW  DATA  PRESS  KEY  *  1" 

4380  PRINT  "" 

4370  PRINT  "            TO  REVIEW  PREVIOUSLY  REDUCED  DATA  PRESS  KEY  I  2" 

4400  PRINT  *   " 

4410  ON  KEY  1  IAPFI  '  NFW  DATA"  COTO  New 

4428  ON  KEY  2  LABEL  '  OLD  DATA"  GOTO  Review 

4430  Wait:  GOTO  Wait 

4445  New:  PRINT  CHR$(12) 

4450  PRINT  "        PUT  THE  DISK  IN  THE  LEFT  DRIVE  IF  IT  IS  NOT  ALREADY" 

4460  PRINT  " 

4470  PRINT  "          ENTER  THE  NAMES  OF  THE  FIRST  AND  SECOND  FILES." 

4480  PRINT  " 

4470  PRINT  "           EACH  FILE  HAS  DATA  FROM  BOTH  DIODE  ARRAYS." 

4500  PRINT  "" 

4510  PRINT  "        FIRST  IS  NO-PARTICLES  I  NEXT  IS  PARTICLES  D1$,D2$" 

4520  PRINT  "" 

4530  Old=0 

4548  SUPEXIT 

4550  Review:  01d=l    'THIS  VARIABLE  IS  PASSED  TO  THE  MAIN  PROGRAM  TO  INDICATE 

4560  • THAT  THE  DATA  TO  BE  READ  TR  ONE  (1)  SCAN  AND  THAT 

4570  'THE  AVERAGING  ROUTINES  WILL  NOT  BE  USED 
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4580 

3!  SPEND 

459(1 

SUB  Readata 

4600 

COM  /ReadaW  B,P,H,Q3*[20], 94*1201, Zz$ 

4610 

Xyz=l 

4f20 

P=8 

4630 

R 1=65536 

4640 

IF  P=2  THEN  Xyz=4097 

4650 

IF  P=2  THEN  Rl=32763 

4660 

IF  P=2  THEN  B=4 

4670 

PRINT  "  • 

4680 

PRINT  n                READING  D 

4690 

ASSIGN  PDi^kl  TO  Q3MtZil 

4700 

ASSIGN  BPisk2  TO  Q4*iZz* 

4710 

ASSIGN  PBtiffl  TO  BUFFER  Yl(») 

4720 

ASSIGN  ?Btif*2  TO  BUFFER  Y2U) 

4730 

CONTROL  PDi?kl,5;Xyz 

4740 

TRANSFER  BDiskl  TO  BBuffljCOUNT  Rl 

4750 

WAIT  FOR  EOT  PDiskl 

4760 

4770 
4780 

CONTROL  PDisk2,5;Xyz 

TRANSFER  P0isk2  TO  BBuff2;CQUNT  Rl 

4790 

WAIT  FOR  EOT  PDisk2 

4800 

ASSIGN  BDiskl  TO  R 

4810 

ASSIGN  BDi5k2  TO  * 

482P 

ASSIGN  BBoffl  TO  * 

4830 

ASSIGN  BBuff2  TO  * 

4840 

SUBEND 

4850 

SUB  Plot2(F,D(*),X) 

4868 

DIH  Titlel$l251,Theta$l20] 

4870 

GINIT 

4880 

IF  X=l  THEN  PLOTTER  IS  705,'HPGL' 

4898 

DEG 

4900 

GRAPHICS  ON 

4910 

VIEWPORT  35,125,35,85 

4920 

Max=10*INT((MAX(D<»))+10)/10) 

4930 

WINDOW  O(l,2),D<E,2),0,Max 

4940 

C=(D(E,2)-DU,2>)/(E-1>*4 

4950 

F=(INKE/4)-l)*4 

4960 

IF  F=0  THEN  F=4 

4970 

CLIP  D(4,2)-G,D(F,2)+2»G,0,Max 

4980 

GRID  G,5,D(4,2)-G,0 

4990 

CLIP  OFF 

5000 

LCRG  8 

5010 

LDIR  90 

5028 

CSIZE  4,. 5 

5030 

FOR  1=4  TO  E  STFP  4 

5040 

HOVE  D(I,2),0 

5050 

LABEL  USING  \DDDD";D(I,2) 

5060 

NEXT  I 

i  FILE  POINTFR  VARIABLE 

i  NUMPFR  OF  SCANS  IN  MOTO»  DATA 

•NUMBER  OF  BYTES      OF  MOTOR  DATA 

! RECORD  *  WHERE  4  SCANS  DATA  BEGINS 

! NUMBER  OF  BYTES      OF  4  SCAN  DATA 

•NUMBER  OF  SCANS 


•OPEN  I/O  PATHS 


iRFADS  NO-PARTICIE  DATA 

'  NOT  AN  OVERLAPPING   TRANSFER 

iXyz  SETS  DISK  FILE  POINTER  TO 

(EITHER  MOTOR  OR  EXHAUST  DATA 

! READS  PARTICLE  DATA 


!I/0  PATHS 


•DEGREES  for  LABEL  DIRECTION 


!AN  X  GRID  LINE  EVERY  4th  POINT 


'MAKES  GRID  UNIFORM 


•PUTS  NUMBERS  ON  X  AXIS 
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11 

LDRG  9 

5090 

FOR  1=10  TO  Max  STEP  10 

5100 

MOVE  D(4,2)-G,I 

5110 

LABEL  USING  ■*,K,jI 

5120 

NEXT  I 

5130 

CSIZE  6,, 6 

5140 

Titlp$='TUO-ANGLE  METHOD" 

5150 

Titlel$="For  Various  Anqle  Ratios" 

5160 

Size$='SIZE  (nicnns)' 

5170 

Thta*='THETA  (rad)" 

5180 

S«b*="r 

5190 

LOIR  90 

5200 

LPRG  5 

5210 

B=D(4,2K-(D(E,2)-D(1,2))/10 

5220 

MOVE  B,H«/2 

5230 

LABEL  Size* 

524  J 

LPIH  8 

5250 

A=(D(E,2)+D(l,2»/2 

5260 

MOVE  A,-Na</4 

5270 

LABEL  Thtal 

528i 

LORG  3 

5290 

MOVE  A,-Hax/4 

5300 

LABEL  S»b$ 

5310 

L3RG  5 

5320 

MOVE  A,H3X*1.15 

5330 

LABEL  Titles 

5340 

CSIZE  4.5 

5350 

MOVE  A,Max«1.05 

5361 

LABEL  Title!* 

5370 

FENUP 

538  J 

5UBEMD 

5390 

SUB  Distribution^**), Tratio,E) 

5400 

LINE  TYPE  1 

5410 

CLIP  ON 

5420 

Tratiot=VAL$(Tratio) 

5430 

MOVE  0(1, 2), D(1,l) 

5440 

FOR  1=1  TO  E 

5450 

DRAW  D(I,2),D(I,1) 

5460 

NEXT  I 

5470 

LHRG  4 

5430 

LINE  TYPE  1 

5*90 

CSIZE  4,. 3 

5500 

MOVE  D(I-l,2),D(I-l,l)+2 

5510 

LABEL  Tratio* 

5520 

PENUP 

5539 

SUBEND 

5540  SUB  Twoangle 

i NUMBER  Y  AXIS 


! STRINGS  FOR  PLOTS 


!  ANGLE  RATIO 

!THIS  SUBROUTINE  PLOTS  THF 
i PART  I  CLE  SIZES  PER  I VED  USING 
'VARIOUS  ANGLE  RATIOS  APPLIED 
!TO  THE  DATA  OVER  A  RANGE  OF 
! ANGLES 
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5550  |xxxxxKX*x»xx*xm****m»»mxxxx#*»*»»*x»»*xxxxx*xx*x»xxj<:»»*x**x 
5560  'x*xx»xxxmm»    SUPPPOGRAM  'TUUANGLE'    **«***#*****»»*»#* 

5570  ixxxxx-xx********.  PARTICLE  SIZING  BY  xkxx»x*k»x*xxxxx*xk 

5590  >**s*««i»x««xx*»  MEASURING  THE  RATIO  k##m***hk««M** 

5590  ixxxxxxxxxxxxxx*  OF  INTENSITY  AT  xx»xx»»»x»x«xxxxxx* 

5600  i*x»xxx*x*xmxx  TWO  ANGLES  X*SW***«KX«#*IK* 

r.h\  0   I  *<*  XXXXXXXX  *»  X»X»«X«*Xft*«ftX«»««»«ftX  X*»»»»«****»«»»»»»K*»**K»»X  XX  X 

5620  OPTION  BASE  I 

5630  COM  /Two/  Avid  024), M,m,F 

5640  COM  /Gauss/  T1(1024),L  'THE  GAUSS  IS  NOT  USED  HERE 

5650  !PUr  THE  CM  PLOCK  HAS  THE!  A 

5660  DIM  D(200,2)  IAND  WAVELENGTH 

5670  PRINT  CUR « (12) 

5680  PRINT  USING  '///////' 

5690  PRINT  '    THIS  SUBPROGRAM  USES.  THE  TWO-ANGLE  METHOD  DESCRIBF.D  BY  BUCHFIE" 

5700  PRINT  "  ■ 

5710  PRINT  "    TO  CALCULATE  PATICLE  SIZE  FOR  VARIOUS  ANGLE  RATIOS  AMD  ANGLES.1 

5720  PRINT  ■  " 

5730  PRINT  "        IT  IS  HOPED  THAT  THE  CURVES  WHICH  RESULT  WILL  SHOW  ■ 

5740  PRINT  "  ■ 

5750  PRINT  ■  WHICH  SIZE  IS  THE  MOST  PROBABLE." 

5760  PRINT  '  ■ 

5770  PRINT  "AFTER  NOTING  FROM  THE  RAW  DATA  WHICH  ANGLES  CONTAIN  THE  CFNTER  LOBF" 

5781  FRINT  ■" 

5790  PRINT  "     ENTfR»*«»*»»*****»THE  SMALLEST  USEABLE  ANGLE," 

5900  PRINT  "  THE  SMALLEST  ANGLE  RATIO, AND" 

5810  PRINT  ■  THE  STFP  SIZE  BETWEEN  ANGLE  RAIOS" 

5320  PRINT  '  YOU  WISH  TO  EXPLORE" 

5830  PRINT  "    EXAMPLE*****»*»m#*      .012, 1.2,. 4" 

5840  INPUT  "  ENTER         thetal ,theta-ratio,theta-step",g,A;B 

5850    OFF  KEY 

5860    X=0  ! GRAPH  WILL  APPEAR  ON  SCREEN 

5970  Begin;  IRUNNING  CONTINUES  HERE  WHEN  A  HARD  COPY  IS  DESIRED 

5880       !Van  Dp  H»lst  and  Gvnprech  &  Sleepevich  exolain  that  the  chanoe  in 

5890       'wavelength  of  the  bea«  is  accounted  for  by  dividing  by  the  index 

59^0       'of  refraction  of  the  nediun.  »***i*«*#L=L/Hl  x»xxxxxx*» 

5910  C=(L/Ml/.57/PI)'2       ! see  page  15  of  nasa  terh  o^oer  2156 

5920  FOR  N=l  TO  1024  !to  see  this  is  a  convenient  constant 

5930     IF  Qai(N)  THEN  5950  IFINDS  POSITION  OF  MINIMUN  ANLGE 

5940  NEXT  N 

5950  FOR  Tratio=A  TO  3  STEP  B  ! VARIOUS  ANGLE  RATIOS 

5960    FOR  J=N  TO  1010/Tratio  STEP  10  'SETS  THE  RANGE  OF  POSSIBLE 

5970  !  ANGLES  (THETA  1) 

5990       Thtal=TKJ) 

5990       Thta2=TKJ*Tratio)  !  THETA  2 

6000       Dpltatheta=Thta2*2-ThtaP2 

6010       11=0 

6020       12=0 
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i0-        ™  lWMh  'mAMWMM m 

6058  NEXT  I  !OF  A  STEADIER  CALCULATION 

6060  11=11/11 

6070  FOR  I=INT(J»Tratio-5)  TO  INT(J»Tratio+5) 

6080  I2=I2+Avl(I)  !  DONE  HERE  FOR  THETA2  FOR 

6C90  NEXT  I  !  THE  GIVEN  ANGLE  RATIO 

6100  12=12/11 

6110  E=(J-N)/10+1  'THIS  IS  A  COUNTER  FOR  THF  AR9AY  CONTAINING 

6120  "PARTICLE  SIZE  AND  THETA1  FOR  THE  GIVEN 

6130  'ANGLE  RATIO 

6140  IF  IKI2  OR  11=0  OR  I2<=0  THEN  D(E,1)=0  'ALLOWANCE  FOR  IF  THE  DATA 

6150  IF  11(12  OR  11=0  OR  I2<=0  THFN  GOTO  Xconp     IIS  NOT  WELL  BEHAVED 

6160  D(E,l)=S9R(-C/Deltatheta«LOG(I?/I1))  IU'0  ANGLE  METHUD 

6170  (VALUE  OF  DIAMETER  based  on  INTENSITY  RATIO 

6190  'For  a  given  ANGLE  RATIO  and  ANGLE  THETA1 

6190  Xconp:      D(E,2Whtal 

6200  Spar=PI»D(E,l)*H1/L         ITHIS  IS  PI»D/LAMRr>A    -  THE  SIZE  PARAMETER 

6210  Tbar=Spar*Tbta2  "THETA  PAR  FOR  THE  LARGER  ANGLE 

6220  IF  Tbar>4  THEN  J=1010/Tratio  ITHIS  ENDS  THE  DO  LOOP  FOR  THIS 

6230  IF  E=l  AND  Tbar)3  THEN  6330  ! ANGLE  RATIO  SINCE  THE  GAUSSIAN 

6240     ! the  above  line  ends  all  calculation  if  IIS  NOT  VALID  WHEN  Tbar  >  3 

6250     'the  first  elewnt  (E=l)failed  the  test 

6260         NEXT  J 

6270  !   FIRST  TIME  THROUGH--  D  HAS  THF  HOST  ELEMENTS  IT  WILL  HAVE 

6290  '   SO  SET  UP  THE  GRAPH  USING  D's  PRESENT  PARAMETERS 

6290  PRINT  CHR$(12) 

6300    IF  Tratio--A  THrN  CALL  Plot2(E,D(»),X) 

6310    CALL  Distribution(D(»),Tratio,E) 

6320  NEXT  Tratio 

6330    ON  KEY  2  LAPEL  "HARD  COPY"  GOTO  Hard 

6340    ON  KEY  3  LABEL  "HENUE  '  GOTO  Subexit 

6350  PRINT  USING  "/////" 

6368  PRINT  '  YOU  CAN  GET  A  HARD  COPY  BY  PRFSSING  KEY  *  2" 

6370  PRINT  "  OR  EXIT  THIS  ROUTINE  BY  PRESSING  KEY  *  3" 

6380  StandbviGOTO  Standby 

6390  Hard:X=l  !A  VARIABLE  TO  CONTROL  PLOTTER 

6400     GOTO  Pea  in 

6410  Subexit:  GINIT 

6420       GRAPHICS  OFF 

6430  SUBEND 

6440  SUB  Shift 

6450  COM  /Readata/  B,P,H,Q3S[20],Q4*l201,Zz$l201,Yl(«)  BUFFER, Y2(»)  BUFFER 

6460  DIM  E(1:BI92) 

6470  PRINT  CHRt(12) 

6480  PRINT  USIKG  "///////" 

6498  PRINT  ■  RAU  DATA  IS  BEING  SHIFTED  TO  CORRECT  FOR  SMALL  GAPS  BETWEEN  SCANS" 

6500  PRINT  " 

6510  PRINT  "  BE  PATIENT,  IT'S  A  LONG  SET  OF  DO-LOOPS' 
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6520 
6530 

6540 

6550 

6561 

6571 

6580 

6590 

6600 

6610 

6620 

6639 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6728 

6730 

6740 

6750 

6768 

6770 

6780 

6790 

6800 

6810 

6920 

6830 

6840 

6850 

6860 

6870 

6880 

6890 

6900 

6910 

6920 

6930 

6940 

6950 

6960 

6978 

6980 

6990 

7009 


THERE  ARE  SOKE  SMALL  GAPS  BETWEEN 
SHIETS  THE  DATA  SO  THAT  THE  FIRST 


TO  THE  VERY  BEGINNING  OF  ITS  1024  FLOCK  IN  THE 


! EXHAUST  DATA  HAS  4  SCANS 


"MOTOR  DATA  HAS  8  SCANS 


D 
OVERALL  ARRAY 
.THE  FIRST  SET  IS  RIGHT  ON,  THE  NEXT  IS  ONE  OFF.THE  THIRD  IS 
TWO  OFF  ,,SO  FORTH.  THIS  MAY  NOT  MATTER  UTTH  OUR  RESOLUTION 
AND  IS  PROBABLY  DUE  TO  THE  MEMORY  CARD  CYCLING  AT  THE  END  OF 
SOME  SCANS.  SEE  THE  CIRCUIT  TIMING  DIAGRAM  IN  THE  THESIS. 
SELECT  B 
CASE  4 

H=3 
CASE  8 

11=7 

END  SELECT 

FOR  K=l  TO  2 

IF  K=l  THEN  MAT  E=  Yl 

IF  K=2  THEN  MAT  E=  Y2 

FOR  J=0  TO  M 

IF  M=7  AND  J<=3  THFN  6800  "One  4096  Block  Doesn't  Need  Shiftino 
FOR  I=(J)»1024+1  TO  (J+1X1024         BLOCKS  OF  1024 
IF  M=3  THEN  L=I+J     lArray  B  Has  the  worst 

•Problem  With  Shifting  Data 
IF  M=7  THEN  L=I+1     iArrav  D  is  Allways  off  bv  one 
IF  LHJ+1)*1024  THENL=M*1)*1034 
I  JUST  TO  AVOID  PROGRAMMING  ERROR  AT  THE  END  OF  THE  APRAY 
E(I)=E(L)  i THIS  SHIFTS  THE  DATA 

(DEPENDING  ON  'I  ', ARRIVED  AT 
!BY  OBSERVING  RAW  DATA 


IOPPFRATES  ON  NO-PARTICLE  AND  PARTICLE  SETS 


NEXT  I 
NEXT  J 
IE  KM  THEN  MAT  Yl  = 
IF  K=2  THEN  MAT  Y2= 
NEXT  K 
EUBEND 
SUB  Store 
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THF  LARGE  ARRAY  CONTAINING  MULTIPLE  SflANS  HAS  BEEN  REDUCFD 
i  TO  A  MEAN  SCATTERING  PROFILE  BY  AVERAGING  AND  FILTERING,  IF 
!  YOU  FEEL  THAT  STORING  THIS  DATA  IS  NECESSARY  THIS  ROUTINE 
!  DOES  IT.  IT  SAVES  DISK  SPACE  TO  STORE  THE  RESULTS  THEN 
!  EUMINATE  THE  RAW  DATA  IF  YOU  FEEL  CONFIDFNT  THAT  THE 
'  REDUCTION  IS  THE  BEST  THAT  CAN  BE.  IN  OTHER  WCRDS, 
i  DO  NOT  PURGE  A  RAW  DATA  FILE  UNLESS  YOU  ARE  ABSOLUTELY  SURE 
'  YOU  WON'T  WANT  TO  REDUCE  IT  AGAIN. 

COM  /Hrdqaits/  Av2<»)  !Av2(»)  IS  THE  REDUCED  DATA 

DIM  Data* 1:1024)  BUFFER  ! BUFFER  USED  TO  TRANSFER  TO  DISK 

MAT  Data=  A«2 

PRINT  CHR$(12) 

PRINT  USING  ■//////' 

PRINT  ■   A  SUGGESTED  METHOD  OF  NAMING  REDUCED  DATA  FILES  IS  AS  FOLLOWS" 

PRINT  '  M MOTOR  BEAM" 

PRINT  ■  X EXHAUST  BEAM" 
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7(110  PRINT  "  C CALIBRATION  ,  IF  NO  'C  THEN  AN  ACTUAL  FIRING  IS  ASSUMFD" 

7020  PRINT  "  G GLASS  BEADS  " 

7030  PRINT  ■  IF  NO  'C  THEN  'G'  STANDS  FOR  'GAP'  PROPFLLANT' 

7040  PRINT  ■  A ALUMINUM  OXIDE* 

7050  PRINT  ■  PPP CHAMBER  PRESSURE  FOR  RUN  OR  OTHER (NOZZLE  TYPE)" 

7060  PRINT  '       PARTICLE  SIZE  FOR  CALIBRATION' 

7070  PRINT  '  MMDD MONTH, DAY  OF  RUN  OR  CALIBRATION' 

7PB0  PRINT  " 

7090  PRINT  "  EXAMPLE:   MCA125JN12' 

7100  PRINT  " 

7110  PRINT  ■  MOTOR  BEAM  CALIBRATION  USING  ALUMINUM  OXIDE  FROM  1  TO  25  HICRONS" 

7120  PRINT  '  ON  JUNE  12" 

7130  PRINT  '  EXAMPLE:   XG550AU10" 

7140  PRINT  ,H 

7150  PRINT  •  EXHAUST  BEAM  DATA  FOR  GAP  P&QPELLANT  AT  550  psi  ON  AUGUST  10" 

7160  PRINT  '  PLACE  A  DISK  IN  THE  RIGHT  HAND  DRTVE" 

7170  PRINT  '  ENTER  THE  FILENAME  YOU  WISH  TO  USE  FOR  THIS  REDUCED  DATA" 

71S0  INPUT  A$ 

7198  CREATE  BDAT  At, 512,16 

7290  ASSIGN  Misk  TO  A* 

7210  ASSIGN  Muff  TO  BUFFER  Data!*) 

7220  CONTROL  BBoff,3;l, 8192,1  (BUFFER  IS  FULL 

7230  TRANSFER  Muff  TO  Misk;COUNT  8192 

7240  WAIT  FOR  EOT  gDisk 

7250  ASSIGN  Muff  TO  * 

7260  ASSIGN  Misk  TO  » 

7?70  SUPENJ) 

7280  PUB  Re«ieu(Avl(«)) 

7290  DIM  Data( 1 : 1024)  BUFFER 

7300  PRINT  '  EACH  REDUPED  FILE  CONTAINS  ONE  SCAN  OMY.  YOU  MUST  KNOU  IF  IT" 

7310  PRINT  "      IS  EXHAUST  ,  MOTOR  CAVITY  ,  OR  CALIBRATION  DATA." 

7320  PRINT  " 

7330  PRINT  '  THE  DISK  WITH  THE  REDUCED  FILE  SHOULD  BE  IN  THE  RIGHT-HAND  DRIVE" 

7340  PRINT  '" 

7350  PRINT  '  ENTER  THE  FILENAME  OF  THF  REDUCED  DATA" 

7360  PRINT  "  TO  BE  REVIEWED' 

7370  PRINT  " 

7388  TNFUT  A* 

7390  ASSIGN  Misk  TO  A$ 

7400  ASSIGN  Muff  TO  BUFFER  DataU) 

7410  CONTROL  Misk,5;l 

7420  CONTROL  Muff ,3;1 ,0,1  (THIS  IS  AN  EMPTY  BUFFER 

7430  TRANSFER  8Di<k  TO  Muff {COUNT  8192  I1024*8=NUM9ER  OF  BYTES 

7440  WAIT  FOR  EOT  PDisk 

7450  ASSIGN  Misk  TO  * 

7460  ASSIGN  Muff  TO  » 

7470  HAT  A»l=  Data 

7488  SUBEND 

7490  SUB  Filter(A(»),Fil) 
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m 

DIM  £(1:1024) 
PRINT  CHP$(12) 

7520 

PRINT  USING  "////////' 

7530 

PRINT  " 

7540 

FOR  J=l  TO  Fil 

7550 

FOR  1=1  TO  1 

7560 

B=I+1 

7570 

C=I+2 

7580 

D=I+3 

7590 

F=I+4 

7600 

G=I+5 

7610 

H=I+6 

7620 

K=I+7 

7630 

L=I+8 

7640 

N=I+9 

7650 

P=I+1 

7660 

E(G)=(A(I)+A(B>+A(C)+A(D) 

7670 

NEXT  I 

7680 

m  a=  f 

7691 

NEXT  J 

7710 

SUBEND 

FILTERING  THE  SCATTERING  PROFILE" 


!THIS  IS  A  Sywetric  Having  Avpraqe 
•TYPE  OF  DIGITAL  FILTER,  EACH 
•DATA  POINT  IS  EQUALLY  WEIGHTED 
!IN  THIS  CASE  BUT  THIS  CAN  BE 
! CHANGED  IF  ONE  DETERMINES  THAT 
"FEWER  P0TNT3  WITH  UNEQUAL  WEIGHTS 
! WOULD  BE  FASTER  OP  GIVE  BETTER 
•RESULTS.  THIS  TYPE  QF  FILTER  WAS 
!U3ED  SINCE  IT  INTRODUCES  NO  PHASE 
•LAG      (Angular  Resolution) . 
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